Cepheids and Long Period Variables in M 33 
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Abstract. We are conducting a long-term photometric survey of the nearby galaxy M 33 to discover Cepheids, eclipsing 
binaries, and long-period variables. The dataset combines previously-obtained optical images from the DIRECT project with 
new observations acquired at the WIYN 3.5m telescope. The entire data set spans over 7 years with excellent synoptic coverage 
which will enable the discovery and characterization of stars displaying variability over a wide range of timescales (days, 
weeks, months, years). 

In this preliminary work we show representative light curves of different variables we found so far in two fields, color- 
magnitude diagrams, and optical Cepheid Period-Luminosity relations for M33. The ultimate goal of the project is to 
provide an absolute calibration of the Cepheid Period-Luminosity relation, and to study its metallicity dependence at optical 
wavelengths. 
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THE PROJECT 

Cepheid variables are a fundamental distance indicator 
that can be applied out to ^40 Mpc with current instru- 
mentation. Our project aims to provide an absolute cal- 
ibration and a determination of the impact of metallic- 
ity on the Cepheid Period-Luminosity (PL) relation, also 
known as the Leavitt Law. Approximately 900 BVI im- 
ages of the Local Group galaxy M33 were obtained be- 
tween 1996 to 2001 by the DIRECT project, mainly us- 
ing the F. L. Whipple Observatory 1 .2m telescope. An- 
other ~1000 BVI images were obtained in 2002-2003 
with the WIYN 3.5m telescope to extend the time cover- 
age and provide higher spatial resolution. The full dataset 
covers the whole disk and spans about 7.4 years. This 
will enable us to identify Cepheids across the disk of 
M 33, which is known to display a large metallicity gra- 
dient 111, |2[] . With such a long time coverage we also ex- 
pect to discover many other variable stars, including long 
period variables. 

In this contribution, we present a preliminary analysis 
of two fields. Representative light curves of Cepheids 
and long-period variables as well as the color-magnitude 
diagrams and PL relations for those fields are shown. 

DATA AND PHOTOMETRIC ANALYSIS 

We analyzed BVI ground-based images that cover two 
fields in M33. The first field is a 10'x5' field centered 
on the spiral arm located north of the nucleus (hereafter 
mOb field). The second is a 11.5'xll.5' field west of the 



nucleus centered on the western arm (hereafter Y3ZC 
field). The Y3ZC dataset originates from the DIRECT 
projects and spans 800 days. The mOb field dataset in- 
cludes WIYN observations in addition to the DIRECT 
data and therefore covers over 2700 days. 

PSF photometry was performed using DAOPHOT and 
ALLFRAME in the three bands 03J. The photometric 
zeropoints and astrometry were determined using the 
M 33 photometric catalog by Massey et al. Q]. A color- 
magnitude diagram of the combined fields is presented 
in Fig. Q] The variability index, J @1, was calculated to 
identify variable star candidates. We used the CLEAN 
algorithm |6] to determine periods for all stars displaying 
a J-index larger than 0.75. We then fit the observed V- 
and I-band light curves to the templates from Stetson 
I5J to identify Cepheid stars. We then applied the same 
technique to the I-band light curves for stars with a J- 
index larger than 2.0 to find variable stars that were too 
red to produce a reliable light curve in the V-band. 

PRELIMINARY RESULTS 

We found more than 200 Cepheids and over 400 long 
period red variables thanks to the large time interval cov- 
ered by our data. Proportionally speaking, we discov- 
ered about twice as many long period red variables in 
the smaller mOb field due to the extended time cover- 
age provided by the additional WIYN data. Examples of 
Cepheid light curves detected in the Y3ZC field together 
with their best-fit template are presented in Fig. [2] 
In the mOb field alone, we detected about 230 long 
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FIGURE 1. Color-magnitude diagrams for the two fields 
combined. Left: all stars detected in the V- and I-bands. Right: 
variable stars only, with the Cepheids represented by crosses 
and long period variables (> 150 days) by triangles. Long 
period variables with no V-band detection are not included. 
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FIGURE 2. Examples of Cepheid light curves detected in the 
1 1 .5' x 1 1 .5' Y3ZC field, as a function of phase. For each panel 
the B-, V-, and I-bands are the lower, center, and upper curves 
respectively. The circles are the observations and the lines are 
the best-fit templates. Periods are in days. 



period variables (P>100 days). A sample of their typical 
I-band light curves are shown in Figs. [3] and [4] V- and B- 
band light curves were not found for many long period 
variables due to their extremely red colors. Most of them 
display an I-band luminosity variation of at least 1 mag, 
sometimes reaching 2.5-3.0mag for the most extreme 
cases. Because of their location on the color-magnitude 
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FIGURE 3. Examples of I-band light curves from long pe- 
riod variables detected in the mOb field, as a function of phase. 
Periods are in days. 
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FIGURE 4. Examples of I-band light curves from ultra long 
period variables detected in the mOb field, as a function of 
phase. Periods are in days. 



diagram (Fig. [TJ, most of the long period variables are 
likely to be late-type stars. A significant fraction of them 
are seen on the red supergiant branch. The location of 
the long period variables on the Wesenheit dust-free 
PL relation (Fig. [5]) also shows that they are in very 
good agreement with the carbon-rich and oxygen-rich 
PL relations for Miras and other semi-regular variables 
J3l, which is consistent with their large variability of 
1-3 mag. In the future analysis of the whole dataset we 
expect to find many more long period red variables and 
study their PL relations in greater detail. 

We generated the Cepheid PL relations for both fields 
(Figs. [5] and |6]l and fit them with the PL relation from 
the OGLE survey for the LMC (§]. The deviations be- 
low 8 days are most likely due to blended objects and 
overtone plusators. Our preliminary distance modulus for 
M33 is consistent with Macri [9]. However, the analysis 
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FIGURE 5. Wesenheit I index of variables detected in the I-band in both mOb and Y3ZC fields. The lines indicate locations 
of different types of variables as reported by Soszynski et al. |7|) for the LMC. Co and C' Q : Oxygen-rich Miras and semi-regular 
PL-relations. Cc and C' c : Carbon-rich Miras and semi-regular PL-relations. D^: Oxygen-rich long secondary period PL-relation. 
Ceph: Cepheid PL-relation. 
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FIGURE 6. Period-luminosity relation of Cepheids found in 
the two analyzed fields of M33. Dashed lines: lo" deviation. 



of the whole dataset is needed to establish a more precise 
distance to M 33. 



FUTURE WORK 

The ultimate goal of this project is to study the PL 
relation dependence on metallicity by using Cepheids 
detected at various radii. The entire dataset is currently 
being analyzed in a similar way in order to gather a large 
sample of Cepheids and other long-period red variables. 
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